Estimates of rates for dissociative recombination of NO^ + e via various mechanisms 
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We estimate rates for the dissociative recombination (DR) of NO J + e~ . Although accurate 
excited state potential energy curves for the excited states of the neutral are not available, we 
estimate that the 1 2 $ 9 and the 1 2 II 9 states of the neutral may intersect the ground state cation 
potential energy surface near its equilibrium geometry. Using fixed nuclei scattering calculations 
we estimate the rate for direct DR via these states and find it to be significant. We also perform 
approximate calculations of DR triggered by the indirect mechanism, which suggest that the indirect 
DR rate for NO2 is insignificant compared to the direct rate. 
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INTRODUCTION 

NO2 is a fundamental molecule in physical chemistry 
that is present in appreciable quantities in the Earth's 
atomsphcre. It has remained one of very few simple tri- 
atomics for which the dissociative recombination process 
has received virtually no attention from cither experi- 
ment or theory. The present work is intended to pro- 
vide a pilot study of this molecule, to build on other re- 
cent progress in understanding triatomic and polyatomic 
DR Vital information about NO2 has been gleaned 
from extensive pioneering studies by Grant and cowork- 
ers These publications have studied in depth the 
Rydberg spectroscopy of low lying vibrational states and 
found that both predissociation and autoionization pro- 
cesses contribute to the lifetime of these states, in varying 
ratios depending upon the vibrational state in question. 
One key finding in particular [i| is that bending excita- 
tion enhances predissociation. Moreover, these studies 
provide evidence that the Born-Oppcnheimer potential 
energy surfaces of certain valence electronic states of the 
neutral cross that of the ground state of the cation near 
the Frank- Condon region of the latter. 

The characteristics of molecular Rydberg states are 
pertinent to the process of dissociative recombination 
(DR) [Irfzi AB+ + e~ -> A + B. Predissociation of 
Rydberg states by valence states indicates that the di- 
rect mechanism ja, 0| can contribute to DR; the direct 
mechanism is competitive only when some excited Born- 
Oppcnheimer electronic states of the neutral are isoen- 
crgctic with the ground state of the cation within the 
Franck- Condon region of the latter. Such crossings of 
potential energy surfaces are responsible both for predis- 
sociation of Rydberg states and dissociative recombina- 
tion (DR) via the direct mechanism. 

The observation of vibrational autoionization in NO2 
indicates that the indirect mechanism [Icl ] may also be 
important. Most studies thus far on DR via the indi- 
rect mechanism (for example, Refs. [lll - fl^ ) have con- 
centrated on molecules with light atomic constituents. 
We anticipate that that indirect DR proceeds violently 



and may lead to the severing of multiple chemical bonds 
in the target [H, . NO2 is both made of heavy atomic 
constituents and they are quite strongly bonded, having 
double bonds between the nitrogen and each oxygen, and 
is therefore a test case for this picture. 

Recently considerable theoretical progress has been 
made in understanding the indirect mechanism and ac- 



curately calculating rates for it [14|, [17| . We apply recent 



theories |18( to DR of NO^ in reduced dimensionality in 
an effort to quantify any possible contribution of the in- 
direct mechanism to the DR rate. Grant's studies show 
that vibrational autoionization is significant for low lying 
vibrational states of the cation; thus the implication of 
Grant's experiments to the question of DR is unclear. 



STRUCTURE CALCULATIONS FOR THE 
DIRECT MECHANISM 

Several authors have looked at the high lying electronic 
states of NO2, although none has explicitly addressed the 
difficult question of the energy of the high lying neutral 
excited states relative to that of the cation. Katagiri and 
Kato [l9[ performed structure calculations on NO2 in 
the course of a treatment of the photodissociation of this 
molecule. They analysed the six lowest doublet electronic 
states of the neutral (lowest at the equilibrium geometry 
of the neutral). Kurkal et al. [IB] calculated potential en- 
ergy surfaces for the lowest two. Takeshita and Shida [2l[ 
calculated many low-lying electronic states of the cation 
in the course of a study of photoionization of NO2, and 
they also analyzed the vibrational structure. 

The lowest two states are the ground 1 2 A\ state, which 
has a bent equilibrium geometry, and the 1 2 B\ state; 
these both become the X 2 n„ state in &ooh geometry. 
The next four states are relevant to the present study. 
They are the first and second states of 2 B2 and 2 A2 sym- 
metry. The 1 2 A 2 and 2 B 2 become the 1 2 $ g state at 
dooff geometry, and the 2 2 A 2 and 2 B 2 become the 1 2 n g 
state at dooH geometry. 

The ground state configuration of NO^ (using labels 
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FIG. 1: (Color online) Structure calculations for the direct 
mechanism as described in text. 



appropriate to c 2 y point group symmetry) is (la? 2a? 
U2 2 3a? 2b 2 2 4a? 3&1 5a? lb\ 4o 2 la|), just like the C0 2 
molecule; the neutral has the 23rd electron in the 6ai or- 
bital. The 1 2 A 2 and 2 B 2 states are described by single 
excitations from the neutral, la 2 and 4&i — > 6a\, respec- 
tively, while the 2 2 A 2 and 2 B 2 are described by 4&i and 
la 2 — > 2bi, respectively. Relative to the cation, therefore, 
all of these states are described by double excitations 
from the NO^ + continuum, and we expect them to 
exist as Feshbach resonances for regions of nuclear geom- 
etry in which they are electronically unbound. As such, 
they provide good candidates for direct DR. 



STRUCTURE CALCULATIONS 

We have performed a set of bound state electronic 
structure calculations using the COLUMBUS quantum 



chemistry suite [22H26| to analyze the crossings of the 
potential energy surfaces of these states with that of 
the cation. These arc performed at the configuration- 
interaction with all singles and doubles (CISD) level of 
theory and Dunning's aug-cc-pvtz basis set |27j, using 
a complete active space self-consistent field (CASSCF) 
reference space with orbitals obtained from a multicon- 
figuration self-consistent field (MCSCF) calculation. We 
include fourteen valence molecular orbitals, leaving out 
only the highest energy 5o 2 (6er„) orbital; we exclude the 
l-3ai and l-2o 2 orbitals from the active space, keeping 
them doubly occupied, giving CAS(12/9) and (13/9) for 
the cation and neutral calculations. 

For the MCSCF calculation that defines the orbitals 
we minimize a weighted average of the ground electronic 
state cation energy and the energies of the first seven 
states of the neutral (the states listed above plus the 
2 S+). In this calculation, the weighted average of the en- 
ergies of the included states is minimized with respect to 
variations of the of orbital and configuration coefficients. 
The balancing of the (N-l)- and N-electron energies is a 
difficult problem, but crucial to the question of DR via 
the direct mechanism. In order to provide an estimate 
of the uncertainty in these relative energies, we vary the 
weighted average; for the seven neutral states we use a 
factor of 1, and for the cation state we use factors of 3.5, 
7, and 14 in the weighted average. 

The results are shown in Fig. [TJ We find that the po- 
tential energy surfaces of the four high lying electronic 
states listed above - the 1 and 2 2 A 2 and 2 B 2 states, 
which correspond to the 2 $ 9 and 2 n s states at linear 
geometry - appear to have crossings with the cation po- 
tential energy surfaces within the Franck-Condon region 
of the latter. 

This bound state structure calculation produces 
not only eigenvalues approximating the location of 
metastable states but also discretized continuum states 
that are unconvcrged artifacts 28 1. One spurious root 



is removed from our plot of bending potential curves in 
Fig.ffl 



SCATTERING CALCULATIONS FOR THE 
DIRECT MECHANISM 

Scattering calculations a re p erformed using the poly- 
atomic UK R-matrix code [if. The GTOBAS code 
is used to obtain a set of diffuse Gaussians used for the 
exterior part of the wavefunction. The same 14 orbital 
valence reference space defines the target states of the 
cation. These orbitals emerge from an improved vir- 



tual orbital calculation [3l| using neutral Hartree-Fock 



orbitals. The Fock operator of which the improved vir- 
tual orbitals are eigenfunctions is defined as an average 
of that for the neutral and the cation. 

All four states corresponding to the 1 2 n 9 and 1 2< 3> g 
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FIG. 3: (Color online) Direct DR cross section versus incident 
electron energy calculated for DR via the 1 



$ 9 state. 



FIG. 2: Location (top) and width (bottom) of the 1 2 B 2 (1 
2 $ g ) Feshbach resonance for linear geometry as a function of 
bond lengths. 



states at linear geometry are observed as Feshbach res- 
onances in these calculations. As might be expected, 
because they couple to the p-wavc continuum, the com- 
ponents correlating with 2 i?2 point group take on a larger 
width than do their 2 A2 counterparts, which couple only 
to partial waves with I = 2 and higher in C2V symmetry. 
As discussed above, the lower energy 1 2 $ 9 state (the 
1 2 A2 and 1 2 B2 states) becomes isoenergetic with the 
cation closer to the equilibrium geometry of the neutral 
than is the 1 2 II„ state. 



ESTIMATE OF DIRECT DR CROSS SECTION 
USING THE REFLECTION PRINCIPLE 

It is possible estimate the magnitude of the direct DR 
rate for NO2 + e~ via the components of the 1 2 $ g state, 
which crosses the cation state closer to its Franck-Condon 
region. To accomplish this, three-dimensional complex- 
valued splined representations of the location and width 
are developed for this 2 B2 Feshbach resonance. The lo- 
cation and width of the calculated resonance at linear 
geometry are shown in Fig. [2] 

The present estimate is based on the standard first or- 
der perturbative expression Q for the DR cross section, 
with the assumption that the survival probability is close 
to unity here. We use the reflection principle (32| approx- 



imation to the bound-continuum matrix element, 

1 f 



dR S [E - (T) - V(R) ^ l {R)V- a 



(1) 



(2) 



The ground vibrational state is approximated as a 
Gaussian in the normal modes but the asymmetric 
stretch has been omitted despite having calculated three 
dimensional surfaces. This omission allows us to treat the 
2 A2 and 2 B2 components of the 1 2 $ s state separately. 

Fig. Q] shows results using a two dimensional vibra- 
tional treatment in both the bend and symmetric stretch 
coordinates, and a one-dimensional treatment including 
the bend degree of freedom only. The static degree(s) of 
freedom are fixed at their equilibrium values. Evidently, 
the cross sections via 2 B2 are higher; this follows from 
the larger electronic decay width of this state. 

We expect the bending degree of freedom to be most 
important for the direct DR mechanism, because it is in 
that degree of freedom that we find crossings of the neu- 
tral and cation states close to the minimum of energy of 
the Born Oppcnhcimcr potential energy surface. Indeed, 
the cross sections obtained from the treatment with only 
the bending mode are similar to those calculated using 
two degrees of freedom, validating this hypothesis. These 
results should be regarded as approximate because there 
is considerable uncertainty in the exact location of the 
cation-neutral crossings. 



INDIRECT MECHANISM FORMULATION OF 
DIRECT DR VIA RYDBERGS 

We follow the method originally outlined in Ref. fljj |. 
This method uses Sicgcrt states in place of the bound 
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FIG. 4: (Color online) Thermally averaged rate coefficient for 
direct DR via the 1 2 $ 9 state. 

vibrational states of the usual vibrational frame transfor- 
mation. These vibrational states Xi are n °t orthogonal 
with respect to the usual hermitian inner product. In 
the present case, as in our previous work on the high en- 
ergy peak of HeH + [l7| , we include cation channels that 
are energetically open. An improved formula is used to 
calculate the unitary defect Di of the S-matrix for an 
electron incident in channel i , not 

°i =i -ei^- r ( 3 ) 

3 

but 

D, I ^.SJ',^ (4) 

jk 

U jk = J dR X i(R)Xj(Ry ■ (5) 

Following the previous work using our method we define 
the DR cross section as proportional to the unitary defect 
of the frame-transformed S-matrix. 

Molecular orbitals are calculated for the cation using 
the MOLPRO quantum chemistry suite j3^|. The cal- 
culation treats 16 electrons in 11 orbitals, with only the 
Is orbitals constrained to be doubly occupied, and min- 
imize the average energy of the 1 A', 3 A', and 3 A" states 
with respect to variations in the orbital and configuration 
coefficients. The first is the ground state at the equilib- 
rium geometry. The next are the lowest triplet states A' 
and A" symmetry Within the Franck-Condon region in 
dooff geometry, these are the 1 3 £~ and one component 
of the 1 3 A U state. These states change discontinuously 
where they are crossed by the 1 3 II g state, each crossing 
being an actual intersection due to symmetry. The po- 
tential energy curves of these states are plotted in Fig [5] 
We use these orbitals in another calculation using the 
UK R-matrix codes, including these three target states in 
the scattering calculation. Wc obtain the i?-matrix from 



FIG. 5: (Color online) Born oppenheimer electronic states of 
NO 2 included in the present indirect calculation: energies for 
r-2=2.12ao, linear geometry, with respect to n. 
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FIG. 6: Total indirect (direct via Rydbergs) probability from 
ID frame transformation calculation using energy- and R- 
independent quantum defects. 



these calculations and solve for the MQDT S-matrix. We 
find that the calculation is very well converged with the 
three target states and a sizable virtual orbital basis in- 
cluded for penetrating configurations. We include s, p, d, 
and /-wave scattering orbitals and obtain a multichannel 
S-matrix with all such partial waves for each of the three 
included target states. 

This treatment is based on the three electronic state 
MQDT S-matrix calculated at zero energy in the 1 *A' 
channel (such that the triplet channels are closed) at the 
equilibrium geometry of the cation. A very simple treat- 
ment is adopted in which the i?-dependence of the quan- 
tum defect is neglected (and the energy dependence is 
likewise neglected). Such a treatment was sufficient to 
yield a sensible calculation of the the high energy DR 
peak of HeH + [T3| • We include only the motion of single 
bond stretch at linear geometry. The results are shown 
in figure [5] 
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CONCLUSION 

We have presented the first estimate of dissociative re- 
combination rates of NO \. Several results analyzing the 
contribution of direct and indirect processes to dissocia- 
tive recombination of NO + e~ were presented. These 
results compared the relative contribution of the direct 
and indirect mechanisms that drive dissociative recom- 
bination. These calculations in reduced dimensionality 
indicate that the direct DR mechanism is likely to domi- 
nate the indirect, the former of which is driven by inter- 
sections of neutral Born-Oppenheimer potential energy 
curves with that of the cation. We find that the 1 2< 3? ff 
and the 1 2 H g electronic states are good candidates for 
states that drive the direct mechanism. There is consid- 
erable uncertainty in the present results, both due to the 
difficulty in locating the aforementioned crossings, and 
due to the reduced dimensionality of the calculations. 
However, these results provide guidance for experimental 
studies of dissociative recombination in these molecules, 
and also provide information that is useful in the de- 
scription of the decay channels of Rydbcrg states of the 
molecule, a topic that has received significant interest in 
the literature. 
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